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1 Introduction 



The nonlinear Schrodinger (NLS) equation in the (1 + 1) space-time dimension is one of 
the completely integrable systems, i.e., the soliton equations [HE]- This model equation 
has been extensively studied to describe nonlinear dynamics in a wide range of physics 
from fiber optics [3llll[5] to Bose-Einstein condensation of cold atoms [HI [71 [8]. The initial 
value problem can be solved exactly via the inverse scattering method (ISM) |9l . In 
particular, the reflection-free condition reduces the inverse problem to a set of algebraic 
equations making it possible to obtain the A^-soliton solution in an explicit way. 

One of major developments in the study of the NLS equation is multicomponent ex- 
tensions preserving the integrability. Manakov [TT] studied a system of the coupled NLS 
(cNLS) equations on the basis of the ISM and obtained the soliton solutions. While the 
interaction of vector solitons for the multicomponent focusing NLS equation is elastic in 
the vector sense it was shown that during the two-soliton collision exchanges among com- 
ponents of each soliton may occur for particular choices of the parameter values [HI [12] . 
In [13], the ISM for a matrix generalization of the NLS equation (in general, in a rect- 
angular matrix form) was developed for solving the initial value problem. By assuming 
the reflection-free condition and vanishing boundary conditions (see below) the A^-soliton 
solution was obtained explicitly. It should be noted that, by appropriate identifications of 
the matrix-field elements, the matrix NLS equation reduces to the cNLS equations of the 
Manakov-type [T2l[T3] and remarkably to the spinor-type that is discovered recently [TH[T5] 
in connection with Bose-Einstein condensates with the spin degrees of freedom. Results 
given in |13J for a general matrix-field, such as the A^-soliton solution, conservation laws 
and Hamiltonian structure, are directly applicable to the reduced systems. Thus, a fur- 
ther analysis of the matrix NLS equation is desired to deal with multicomponent nonlinear 
dynamics under different circumstances. 

In this paper, we study a square matrix NLS equation under the nonvanishing bound- 
ary conditions by means of the ISM. The multicomponent system with such boundary 
conditions is regarded as an extension of a basic single-component NLS equation with the 
self-defocusing nonlinearity studied by Zakharov and Shabat [10] and also that with the 
self- focusing nonlinearity investigated by Kawata and Inoue [16]. As compared to the case 
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with the vanishing boundary conditions [13] , the conservation laws and Hamiltonian struc- 
ture are the same, while the Lax pair, the direct and inverse problems, and the A^-soliton 
solution should be reformulated reflecting the boundary values of the matrix-field. 

This paper is organized as follows. In Sec. [21 nonvanishing boundary conditions for the 
matrix NLS equation are introduced. The Lax pair is provided to formulate the auxiliary 
linear system. Then the conservation laws are constructed systematically. In Sec. [3], the 
direct and inverse problems are solved along the ISM procedure and the multi-soliton 
solution is presented. Section H] is devoted to the concluding remarks. 

2 Formulation 



The matrix NLS equation is expressed as 

iQt + Qxx-2eQQ^Q = (£ = ±1), (2.1) 

where Q{x,t) and O are an / x / matrix valued function and the zero matrix, respectively, 
is the Hermitian conjugate of Q, and the subscripts t, x denote the partial derivatives. 
The case e = —1 {e = +1) of (12. ip is often referred to as the self-focusing (-defocusing) 
one. In [13], it was shown that through the ISM the system has an infinite number of 
conservation laws. The initial value problem was solved and the A^-soliton solution was 
obtained under the constraint e = —1 and the vanishing boundary condition: 

Q{x,t)^0 as x^±oo. (2.2) 

Under these conditions each soliton forms the so-called bright soliton with P components. 
Setting the form 

/ Qii ■ ■ ■ Qii \ 



Q 



Qml 





Qml 





(2.3) 



V ■■■ / 

with m < I, one can achieve a rectangular matrix reduction that is compatible with 
the vanishing boundary condition. In particular, the m = 1 case corresponds to the l- 
component Manakov model [12] . Other types of reduction can be obtained in a nontrivial 
way by putting some components equal without breaking consistency of equations [HI [15] . 
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On the other hand, for another integrable constraint e = +1, leading to 

iQt + Qxx - 2QQ^Q = O, 



(2.4) 



the boundary condition should be altered appropriately, which has not been investigated 
so far. Equation (12.41) is a matrix generalization of the NLS equation for a scalar field 
q{x, t) with a self-defocusing nonlinearity, 

iqt + qxx - 2\q\^q = 0, (2.5) 

which possesses dark soliton solutions. For the self-defocusing NLS equation (12. 5p . the 
boundary condition at x ^ ±oo is assumed to be the nonvanishing one, e.g., a constant, 
\q{x)\ Ao, rather than the vanishing one, \q{x)\ — > 0. The ISM procedure was applied 
to the system ( 12. Sp in [10]. The analysis of the NLS equation under the nonvanishing 
boundary conditions was extended to the self- focusing case in [T6] . 

From now on, we concentrate on the analysis of a full-rank / x / square matrix NLS 
equation. We do not include reductions ( 12. 3p in this work, i.e., the vector (Manakov) NLS 
equation falls outside our considerations. Systematic study of such reductions based on 
the symmetry argument is an open issue. In this section, we introduce a square matrix 
type of nonvanishing boundary conditions for eq. ( 12. ip and formulate the Lax equation for 
the ISM. 

2.1 Nonvanishing boundary condition 

We assume that the / x / matrix valued function Q{x, t) satisfies the following nonvanishing 
conditions, 

Q{x,t) Q± as X — »• ±oo, (2.6) 
QiQ± = Q±Qi = Xli, (2.7) 

where Aq is a positive real constant and I denotes the / x / unit matrix. Noting the freedom 
of unitary transformations: 

Q'=l(QV, (2.8) 
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with W, V unitary matrices, we see that if Q is a solution of eq. (12.11) . then Q' is also a 
solution. By this freedom, without loss of generality, we can fix one side of the boundary 
condition as 



= Aoe'(^"-"*)/. 
To avoid a complexity, separate the carrier wave part, 

Q(x,t) =Q(x,t)e'("), 
where the dispersion relation is determined as 

u} = P + 2eXl 



Then the original nonlinear evolution equation ( 12. ip is equivalent to 

iQt + + 2ikQ, + 2e{Xll - QQ^)Q = O. 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



Here and hereafter except the final expression (I3.100p . we drop the hat of Q for a notational 
simplicity. Accordingly, eq. (12. 9p is rewritten as 



Q+ = Ao/. 



(2.13) 



In what follows, we focus on eq. ( 12. 12^ with the boundary conditions ( 12. 6p . ( 12. 7p . and 
2.2 Lax pair 

We introduce the Lax matrices in the following forms. 



U = iX 



' -I 0' 


+ 





Q ' 


I 








(2.14) 



' -21 


' 







21 


+ x 



2ikl 2Q 
2eQ^ -2ikl 



e{XlI - ggt) + 2ikQ 

si-Ql + 2ikQ^) -siXlI - Q^Q) 



(2.15) 



where A is the spectral parameter that is independent of time, A^ = 0. The potential matrix 
Q satisfies the nonvanishing boundary conditions (12. 6p . (12. 7p with eq. (I2.13p . In the ISM, 
we associate a set of linear problems: 



(2.16) 
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where \1/ is a 2/ x / matrix function. We also use the following representations, 



" ^1 " 


, v = 




Vl2 ' 


^2 


_ ^21 


V22 _ 



(2.17) 



where all the entries are / x / matrices. Then, the Lax equation is obtained from the 
compatible condition of eqs. fl2.16p . 

Ut-V^ + UV -VU = 0, (2.18) 

which is equivalent to the matrix NLS equation f l2.12p . 

2.3 Conservation laws 

In this subsection, we recapture a systematic method to construct local conservation laws 
for the matrix NLS equation [13] . The method was originally developed for the scalar field 
case [17]. 

Introduce the quantity. 



r = ^2^r^- 

From eqs. fl2.14l) - fl2.18l) . one can prove 

|tr(Qr)^|tr(K.r + Kn), 

2iAQr = -sQQ^ + gr, + {QVf. 

Note that eq. fl2.20p has a form of conservation law. Expand QT in A as 



(2.19) 



(2.20) 
(2.21) 



^^"5Z(2iA)r^ 

i=i 



(2.22) 



Then the trace of each coefficient Fj is a conserved density and eq. (12.201) represents an 
infinite number of continuity equations. From eq. (12.211) . we recursively obtain 



Fi = -QQ\ 

tTF2=tT{-QQl), 

trFa = tr(-gQt^ + eQQ^QQ^). 
By direct calculation, all elements of Fi are shown to be conserved densities. 

6 



(2.23) 
(2.24) 
(2.25) 



3 Inverse Scattering Method 



In this section, we carry out the ISM procedure for eq. (12.121) based on the Lax pair fl2.14p 
and fl27[5ll . 



3.1 Direct problem 

We consider the eigenvalue problem, 



eQ^ iXI 



(3.1) 



and define the Jost functions and scattering data for them. Here Q plays a role of potentials 
in the eigenvalue problem. In Sec. 13.11 and 13.21 the analysis will be made with fixed time t. 
Under the nonvanishing boundary conditions (12. 6p . U has the asymptotic forms: 



-iAJ g± 
eQi iAJ 



UiX) - UHX) = 
The characteristic roots of U^{X) are twofold, 



as X ^ ±oo. 



(3.2) 



iC, -iC, 

where ( = (A^ — eXqY^^. We introduce 21 x 21 matrix functions T and by 



(3.3) 



nx,c,x) 



-iQix) (A-C)/ 
(A-C)/ ieQ^ix) 



T±(A,C)= lim T(A,C;x). 



(3.4) 
(3.5) 



Here Q{x) is an / x / smooth matrix function that satisfies the same boundary condition 
as eq. (12.61) . 

Q{x) — Q± as X ±oo, (3.6) 



and for all x. 



Q{x)Q^{x) = Q^{x)Q{x) = Xll. 



(3.7) 



Using T''=(A,C), we can diagonalize U'^i^X) as follows. 



U^{X) = -KT^iX, C)(a^ ® /)[T±(A, Or\ 
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where a* (i = x,y,z) is the Pauh matrix and denotes the direct product, 





a"" (g) 


J = 


'I ' 

-I 








use of eqs. fl3.8|). we define matrix Jost functions ^/^^^ , 




eq. 


(j3.1j), whose asymptotic forms are, respectively, given 


by 




^r~r-(A,c) 


' / " 



e-'^^ as 


X 


— > — oo 




^2-~T-(A,C) 


' ' 

I 


e'^" as 


X - 


— oo. 




V'+~T+(A,C) 


' I ' 




e"'^^ as 


X 


+00 




V'2+~T+(A,C) 


' ' 

I 


e'^" as 


X - 


H> +00. 



(3.9) 



(3.10a) 
(3.10b) 
(3.10c) 
(3.10d) 

We note that {'iptji't} well as {i'1,4'2} constitute fundamental systems of solution. 
This is easily proved by using the usual Wronskian defined by the determinant. In fact, 
one can show 

^det[<l>i,$2] =0, (3.11) 

for any two 21 x / matrix solutions $1, $2 of eq. (13.11) . Checking the value at x ±00, we 
have 

det[^^,^t] = {2C{\-C)y. (3.12) 

This indicates the linear independence of -ipf and tpf except the branch points of i.e., 
A = ±v^Ao at which the solutions degenerate. 

If we use a notation i/j^ = [ip^, t/^^], relations (13.101) can be rewritten compactly in the 
following form, 

V'^(A,C;a;) ^ r^(A,C)J(C^) as x ^ ±00, (3.13) 

where 



JiCx) ^ 

Then the scattering matrix S{\, () is defined by 



O e'^^J 



V-(A,C;a;) = V+(A,C;a;)5(A,C), 
' A{X,() B{X,() 



^(A,C) 



(3.14) 

(3.15) 
(3.16) 
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where all the entries A, A, B, and B represented by / x / matrices constitute scattering 
data. 

It is useful to consider another slightly modified eigenvalue problem [18]. Under a 
t r ansf ormat ion , 



the eigenvalue problem (13. ip becomes 



$^ = [/$, U = T-\UT - T^) = -[(a' I + W, 



(3.17) 



(3.18) 



where W = (Wab), a,b = 1, 2 with / x / matrices, 



Wu = ^{Q^Q + Q^Q-2Xll)- 



2C(A-C)^^ ' 



12 



21 



2X1 



2X1 
is 



Q 



(Qgt + ggt _ 2Xll) + (QQt _ ggt 



c 



(QtQ + QtQ _ 2Xll) + (gtg _ gtQ) 



2C 



w^22 = -^(ggugg^-2A^/) 



2C(A-C)^^^- 



As solutions of eq. ( 13.18p . we introduce new Jost matrices </> 
simpler asymptotic forms, 



(3.19a) 
(3.19b) 
(3.19c) 
(3.19d) 



T with 



4>^{XX', x) ^ J{Cx) as X ^ ±oo. 



(3.20) 



These matrix Jost functions are connected with each other through the same scattering 
matrix (13.161) as 

0-(A,C;x) = 0+(A,C;x)S(A,C). (3.21) 

We now examine the properties of the scattering data. For this purpose, define a matrix 
function for two 21 x / matrix functions $i(A, (; x) {i = 1, 2) as 



M[^,,^2] = ^liX*,C;x 



I o 

O -el 



$2(A,C;a;) 



(3.22) 



If $i(A, x)^ $2 (A, C; x) are solutions of eq. (13. ip . one can show that 

d 



dx 



M[<l>i,$2] =0. 



(3.23) 



From the asymptotic forms (13.101) . we obtain the following relations, 

M[i;t,^t] = -eM[i;t,i;t] = 2eC(A - C)/, 

M[ijt,ijt]=0, 

M[^+ ^r]=2£C(A-CM(A,C), 

M[^2+,7A21 = -2C(A-C)^(A,C), 

M[^2+,7Ar] = -2C(A-C)fi(A,C), 

M[^+ ^2"] = 2<(A-C)5(A,C). 



These relations are rewritten into 
which can be shown, e.g., as 



^(A,C) B{\,0 





'10' 




I 



M[v;r, v^r] =^^(A*, c)M[i,t. A+]^(A, c) + ^^(a*, c)M[i,t, i^tm\ c) 

=2£C(A - C)(^^(A*, C)^(A, C) - eB\X\ C)B{\ Q) 
=2eC(A-C)/- 
The relation (13.251) leads to the inversion of eq. (13.151) . 

V'+(A,C;x) = V-(A,C;x)[^(A,C)]-\ 



-eB^i\*,C) A^{\*,C) 



In this system, we have involution relations for the Jost functions 

V'±(A,C;x)=V^(A,-C;x)P±(A,C), 

where 

P±(A,C) = J{Cx) [r±(A,-C)]"'T±(A,C)J(C 
_ 1 r O ieQi 
A + C [ -iQ± O 

From the involution, we have another set of relations for the scattering data: 

5(A,C)= [P+(A,Cr'5(A,-C)P-(A,C), 
^(A,C) = ^- 



X) 



■g+A(A,-c)gl 



(3.24a) 
(3.24b) 
(3.24c) 
(3.24d) 
(3.24e) 
(3.24f) 

(3.25) 



(3.26) 

(3.27) 
(3.28) 

(3.29) 



^QiBiX,-C)Ql- 



(3.30) 

(3.31) 
(3.32) 

(3.33) 
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By use of eqs. ( I3.15P and fl3.27p . we can prove 

detA(A,C) = detfVT i^t] = det A^{X*,C) 

= (detA(A*,C))*. (3.34) 

Additionally, taking the determinant of both sides of eq. (13.321) . we have 

det A{X, C) = det A{\, -() det [Aq ^Q+QLI . (3.35) 

Combining eqs. (1331) and ^M), we find that detA(A,C) oc (det A(A*, -(*))*. As a 
consequence, if (Aj, Q) is the zero of det A, (A*, — Cj ) is also the zero of det A, and the pairs 
(Aj, —Q) and (A*,Cj) are the zeros of det A. For e = +1, furthermore, the self-adjointness 
of the eigenvalue problem (13. ip leads to \j = A* G M and (j = —(* G iM. 

Finally, we make clear the analytical properties of the Jost functions (13.101) in regard 
to complex A. To this end, we prepare a two-sheet Riemann surface for A where ( = 
(A^ — £:Aq)^/^ is single- valued. For e = +1, cuts are made in (— oo, — Aq] and [Aq, oo) (see 
FigH]). Each sheet is characterized such that Im C > (I^^ C < 0) on the upper (lower) 
sheet. On the other hand, for e = — 1, cuts are made in [— iAo, iAo] (see Figl2]). It is required 
that Im Im A > (Im Im A < 0) on the upper (lower) sheet. The Jost functions satisfy 
the following Volterra-type integral equations, 

0±(A, C; x) = J{Cx) (^I + £ dyJ{Cy)-'Wcf>^{y)^ . (3.36) 

Suppose that 

/oo 
\{Wab) ij\dx < oo, (3.37) 
-oo 

for all a, 6 = 1, 2 and i, j = 1, ■ ■ ■ ,1. We may have the Neumann series solution 
4^^{X,C,x)JiCx)-' = J2 / / dy2--- dy^G{y,)---Giy,,), 

J ±oo J ±oo J ±oo 

= T<(»e^±-^(^^)'^^ (3.38) 



where G{y) = J {({y — x))~^ W J {({y — x)) and T<(>) denotes the time-ordered product. 
Examining the convergence of the Neumann series and its derivatives, it is found that 
0^(A, C; a;)e"'^, 0^(A, C; a;)e~"'^ are bounded and analytic in the region where Im > 0, 
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w 
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and 0^(A, C; x)e'^'^, 02 (-^5 C; a;)e~'^^ are bounded and analytic in the region where Im ^ < 0. 
Relations fl3.24cp and (13.24dp show that A{\,Q (^(-^^C)) is analytic in the region where 
Im C > (Im < 0). We also have the asymptotic behaviors of the Jost functions and the 
scattering data as A, C — ^ oo from the asymptotics of W\ 

'-Q-^Q. + 0{l/\\\), ifC^A, 
0(1/|A|), if C^-A, 

^21 = 

'q,q-^ + o{i/\x\), ifC^A, 

0(1/|A|), ifC^-A. 
In calculating Neumann series, the following formulae are useful: 



22 



(3.39a) 

(3.39b) 

(3.39c) 

(3.39d) 



n-in 



Q-\b)Q{a) 



(3.40) 
(3.41) 



for a < b. The proof is as follows. Write the RHS of (13.401) as X{a). Differentiating X{a) 
gives Xa = XQ~^Qa. Solving this equation, we get X(a) = AQ{a). Since X{h) = I, we 
find A = Q^^{b), and we arrive at the desired formula. We get (13.411) just in the same way. 
Consequently, when ~ A, we have the asymptotics 



0±(A,C;x)J(Cx) 



-1 



o 



O Xo^Qix)Qi 



Aq ^Q+Q- 



o 



o 

Aq '^Q+Q- 



5(A,C) = 

On the other hand, when ( ~ —A, 

0±(A,C;x)J(Cx)-^ = / + O(l/|A|), 
5(A,C) = / + 0(1/|A|). 

Furthermore, one can also show 

^r(A,C;^)e^^1^(A,C)]-'-T+(A,C) 
^2-(A,C;x)e-'^^[A(A,C)]-^-T+(A,C) 



+ o{i/\x\), 

0{1/\X\). 



I 
O 



o 
I 



Im C > 0, 
Im C < 0. 



(3.42a) 
(3.42b) 

(3.43a) 
(3.43b) 

(3.44) 
(3.45) 
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Note that, however, when A = ±y^Ao the scattering data may have a singularity of 
the order 0{1/Q which can be seen from eqs. ( I3.24p . while the eigenf unctions are well- 
defined at those branch points in the general situation [2]. It should be remarked that the 
introduction of Q{x) is irrelevant in the analysis of the original eigenvalue problem (13. ip 
and the inverse problem discussed in the subsequent section. 

3.2 Inverse problem 

In this subsection, we derive the Gel'fand-Levitan-Marchenko equations which give the 
solution of the inverse problem, by use of the Jost functions on the complex Riemann 
surface. 

We assume that the Jost functions (13.101) are expressed as 

V±(A, C; x) = T±(A, C)J{Cx) + / /C(x, s)T±(A, C)J{Cs)ds, (3.46) 

J X 



where the kernel matrix is 



JC(x, s) 



/Cii(x, s) /Ci2(x, s) 

/C2l(x, S) IC22ix,s) 



(3.47) 



with ICij{x,s) being / x / matrix functions. Substituting the expression (I3.46P into the 
eigenvalue problem (13.11) . after some calculations, we obtain a linear differential equation 
for the kernel matrix /C(x, s), 

d^lC{x, s)+{a^ (g) I)dslC{x, s)(cr^ O /) 

+ (a^ ® I)}C{x, s){a' ® /) [f/±(A, C) + iA(a^ ® /)] 

- [U{X, C; x) + iX{a' ® /)] /C(x, s) = O, (3.48) 

with boundary conditions, 

2/Ci2(x,x) = g+ -g(x), (3.49a) 

2IC2iix,x)=eiQl-Q^ix)), (3.49b) 

/Cjj(x, s) — s> O as s ^ ±oo. (3.49c) 

This type of a linear system is known as the Marchenko equations and can be uniquely 
solved, which guarantees the existence of the expression (13.461) . 
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First, we concentrate on e = +1 case. Before going on, we mention several facts. 
For the characteristic roots (13.31) . we have introduced two branch cuts on the real axis, 
(—00, — Ao] and [Ao, 00), as shown in Fig. [H We define contour paths C enclosing a region 
in the upper sheet (Im > 0) of the Riemann surface of A and C enclosing a region in the 
lower sheet (Im C < 0) as (see Fig. [1]), 



c 


=T + B, 


C 


= f + B, 


T 


=r+ + r-. 


f 


= r+ + f-. 


B 


=B+ + B-, 


B 


= B-^ + B- 



(3.50a) 
(3.50b) 
(3.50c) 



where the superscripts + and — denote the part of the path F (F) which exists in the upper 
and lower half plane of each sheet, respectively, and the part of the path B [B) which exists 
in the right and left half plane of each sheet, respectively. The radius of F (F) is assumed to 
be large enough for C (C) to enclose all the zeros of det A (det A). As noticing in Sec. 13.11 
if (Aj, Q) is a zero of det A, it holds that A-,- G M and Q G iM. Correspondingly, (Aj, —(j) is 
a zero of det A. Suppose here that det A and det A have zeros, respectively. Along the 



A' 



B- {B- 



Im A 



V, r+ (f +) 



o ^ti^^) 
• • 



Re A 



An 



F- (f-) 



Figure 1: Cuts, bold lines, and integral contours, dotted lines, on the upper (lower) sheet 
of the Riemann surface of A plane for the self-defocusing case, e = +1. The cross denotes 
the zero for det A (det A) on the upper (lower) sheet. 
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branch cut in the upper sheet, one can show that 



Jb c 



dAe 



ax'- 



0, 



(3.51) 
(3.52) 



where z is real and 6{z) denotes the delta function. Using these formulae, we obtain 



i^^dA^T±(A,C)=5(^)(cr^®/). 



(3.53) 



In the lower sheet, integral formulae fl3.5ip . 03.520 and fl3.53p hold with replacing B ^ B, 
Going back to eqs. fl3.46p . we explicitly write as 



^+(A,C;a;) = T+(A,C)e-^'^^ 
^+(A,C;x)=T+(A,C)e'^^' 
We rewrite eq. (13.151) into 



I 

O 



' ' 




I 





+ / ds/C(x,s)T+(A,C)e-'^ 



X 

oo 



o 



o 
I 



(A, C; x) [A{\ =^i+(A, C; x) + C, x)B{\, () [A{\, ()]-\ 
^2"(A, C; x) [A(A, C)]-^ =^tiX, C; x) + ^+(A, C; x)B{\, c) [A(A, C)]-^ 



/ 
O 



Substituting eq. (I3.54ap into eq. (I3.55al) . we have a relation 
V^r(A,C;a:)[A(A,C)]-'-r+(A,C)e-^^^' ^ 

+ |T+(A,C)e^^^ ^ 

POO 

+ / ds/C(x,s)r+(A,C)e'^'^ 

J X 

When we multiply 



O 
I 



4vrC 



(y > x) 



(3.54a) 
(3.54b) 



(3.55a) 
(3.55b) 



B{X,0[A{XX)r. (3.56) 



(3.57) 



on the both sides of eq. fl3.56p . the left hand side of eq. fl3.56p becomes analytic on the 
upper sheet of the Riemann surface (Im ( > 0), with the exception of the points Xj, at 
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which it has simple poles. Here we have assumed that 1/ det A{X, () has isolated simple 
poles {Ai, A2, . . . , Aiv} in the upper sheet. From eq. (13.441) . we can show that as |A| 00, 
the left hand side of the resultant equation behaves like exp[— Im ({y — x)]C(l/|A|). We 
integrate the relation fl3.56p multiplied by eq. (13.571) along the contour B. In the integration 
of the left hand side, the contour can be closed through infinity, i.e., C = B + T, so that 
the integral of the left hand side is equal to the sum of the residues at A = Xj, 



1 



dX- 



C 



^r(A,C;a:)e^^"[^(A,Cr'-T+(A,C) 



/ 

O 



N 



A{X„Q 



^S^O ^^"^^-^-^^(detAy(A,,o; 



N 



i^e'^^^^2+(A„0;x)n 



J' 



(3.58) 



where A is the cofactor matrix of A and Ilj is the residue matrix at A = Xj defined by 

1 



n. 



Res 

^=^j ,C=Cj 



2C 



i?(A,C)KA,C)] 



-1 



(3.59) 



In the second equality of eq. (I3.58p . we have used a relation deduced from eq. (I3.55ap 
at (A, C) = (Aj,Cj)) which means that ipii^jXji ^) is proportional to ip2i^jy(jj^)- The 
integral in the right hand side of eq. (I3.58P is transformed into 



lC{x,y) 



O 
I 



+ Tc{x + y) 



' ' 




I 





POO 

/ dslC{x, s)J='c{s + y) 

J X 



o 
I 



where 



J^,{z) ~ [ dA^r+(A,C)p(A,C), 
p(A,C) =5(A,C)KA,C)]-^ 



where J-'[c{z) = dj^ic{z)/dz, and 





we 


' ' 




I 





1 

47r 
1 

47r 



00 gi^2 



d^— [P(A,0-P(-A,0] 
dee^«Mp(A,0+p(-A,0]- 



(3.60) 

(3.61) 
(3.62) 

(3.63) 

(3.64a) 
(3.64b) 
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Using the integral representation (13.54bp in eq. (13.581) . we rewrite eq. (13.581) into the form 
(I3.63P as follows, 



' ' 




I 





where 



N 

J•l,(^)=5^i^,e^'^^^ 

AT 

Combining eqs. (13.601) and (13.651) . we finally obtain the integral equation: 



Here 



O 
I 



+ T{x + y) 



O 
I 



+ / dslC{x, s)J^{s + y) 



' ' 




' ' 


I 








Following the same procedure using eqs. fl3.54bl) and fl3.55bl) . we obtain 



}C{x,y) 



where 



/ 

O 



+ f{x + y) 



I 
O 



POO 

/ dslC{x, s)f'{s + y) 

J X 



' I ' 




' ' 











T{z) 



An 



dA^T+(A,C)p(A,C), 



p(A,C) =5(A,C)[A(A,C)] 



-1 



The residue matrix at the zero of det A, say, (A, () = (Aj, Q) is defined by 

1 



Res 



2C 



(3.65) 



(3.66a) 
(3.66b) 



{y > x). 

(3.67) 
(3.68) 

(y > x), 
(3.69) 

(3.70) 
(3.71) 

(3.72) 



Integral equations (I3.67P and (13.691) are called the Gel'fand-Levitan-Marchenko equations. 
Solving these equations as to the kernel }C{x,y) for given J^{z) and ^{z), we can obtain 
the potential matrix Q through the relations (I3.49p . 

Let us move on to e = —1 case. Integral contours are depicted in Fig|2l The branch 
cut in the Riemann surface is made along [— iAo,iAo]. Contours playing the same role as 
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in the previous case are named with the same letter with a prime We define a contour 
path C {C') enclosing a region Im ^ > (Im C < 0) as, 



C =T' + B', C' = f ' + 

r'=r'+ + r'-, r' = r'+ + f'~, 

B'=B'+ + B'-. B' = B'+ + B'~ 



(3.73a) 
(3.73b) 
(3.73c) 



where the superscripts + and — denote the part of paths in the upper and lower half 

Im A 



A' 



B'+ (F+) 



B'- [B'- 



r'+ (f'+) 



iAo 



Re A 



B'- {B'- 



-iAn 



f - (r'-) 



Figure 2: A cut, bold line, and integral contours, dotted lines, on the upper (lower) sheet 
of the Riemann surface of A plane for the self-focusing case, e = —1. The cross denotes 
the zero for det A and det A on both sheets. 

plane of each sheet. The contours B'^ and B'^ are along both the real axis and the branch 
cut. The radius of T' (f') is assumed to be large enough for C (C) to enclose all the 
zeros of det A (det A). In contrast to the previous case, the zero of det A is complex and 
makes a pair. Assume that there are = 2M zeros for det A and label them such that 
A2fc = ^2k-i^ ("ik = —Qk-i for = 1, ■ ■ ■ , M. Correspondingly, we have the zeros of det A 
such that {(Ai, -Ci), (Ai, CDr ■ ■ , {^m, -(m), {^*m, CX/)}- The derivation of the integral 
equation goes in parallel to that for the case e = +1. As a consequence, we arrive at just 
the same equations fl3.67p and (13.691) . 
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We remark on the properties of the matrices Uj and Uj. First, the determinant is zero. 
This is proved directly as 

det5(detA)'-i 



det n, oc , 

^ [(det A) 

= 0, 



7^1 (^J'Ci) 



(3.74) 



where we have used det A{Xj, (j) = 0. Similarly, we have det 11^ = 0. Second, we have the 
following relations: 



n, = n] = -n, 



(s = -l). 



These are proved as follows. From eqs. (13.251) . (I3.3ip . and (13.321) . we obtain 

BA-\\,C)=e[iBA-')iy,C)]\ 
BA-\XX) = -eBA-\X,-C). 



(3.75a) 
(3.75b) 

(3.76) 
(3.77) 



For example, we demonstrate Uj = — Ilj for e = +1. Substituting eq. (13.760 into the 
definition of II j (13.721) we have 



n, = Res 



[{BA-'){y,cW 



2C 



Res 



iBA-'){X*,C' 



-n 



Res 

A=A„C=0 
t 



Y^{BA-^){\X) 



J' 



(3.78) 



where in the third equality we have replaced the dummy variables as A* A, C* — » C and 
used A* = Aj, Q = —Q. The rests are obtained similarly. 

3.3 Time dependence of the scattering data 

Next, we consider the time dependence of the scattering data. Under the nonvanishing 
boundary conditions (|2.6p . the asymptotic forms of the Lax matrix V are given by 



2{\-k) 



2iXI g± 
i 2iXI 



eQ 



as X ±oo. 



(3.79) 
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Operating fl3.79p on the asymptotic forms of the Jost functions (13.131) gives 



I 
O 



o 
I 



g-iC- = _ 2i(A -k)[C + Xia' ® /)]T±(A, C) 

e^^^ =2i(A-fc)[C-A(a^®/)]r±(A,C) 

±{t) 



/ 

O 



o 
I 



Then we define the time-dependent Jost functions 'ipj , j = 1, 2, as, 



±W ^ g-2i{A-A;)[C+A(<7-$5/)]t^± 

^ ^---2i{\-k)[C+X(a-<S>I)]t-Kxrp± 



I 

o 



±{t) ^ g2i(A-fc)[C-A(a-®/)]t^± 

^ g2i(A-fc)[C-A{a-®/)]i+iCxj.± 



which obey, respectively. 



±(i) 



dt 

at 



±{t) 



cLS X 



cLS X 



±CXD, 



±oo, 



±(t) 



These relations give 



dt 
dt 



{V + 2i(A -k)[C + A((T^ ® 
{V-2i{\-k)[(-\ia^®I)]}ij, 



Substituting the definitions of the scattering data (13.151) : 

^r(A, C; X, t) =tPt{K C; x, t)A{X, C; t) + tPtiK C; x, t)B{X, C; t), 
C, X, t) =V'i+(A, C; X, t)B{X, C; t) + ^PtiK C; a;, t)A{\, (■ t), 



(3.80a) 
(3.80b) 



(3.81a) 
(3.81b) 

(3.82a) 
(3.82b) 

(3.83a) 
(3.83b) 



(3.84a) 
(3.84b) 



into eqs. (I3.83p . and then taking the limit x ^ oo, we find the time dependence of the 
scattering data as follows: 



A{X,C,t) =A(A,C;0), 

n,(t) =n,(o)e^'^^(^^ 



(3.85a) 
(3.85b) 
(3.85c) 
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Using eqs. (13.851) . we obtain explicit time dependence of J-'i{z,t) and J-'2{z,t), 

:Fi{z, t) =- J d^- [p(A, 0) - p(-A, 0)] - J] in, (0)e^^^^+^^^^(^^-'=)*, 

(3.86a) 

J^2iz, t) =— / dee'«^+^'«(^-^')* [p(A, e; 0) + p(-A, ^; 0)] - V iA,n,(0)e'^^-+^^^^(^^-^)*. 

(3.86b) 

We have similar time dependence for ^i{z,t) and ^2{z,t) with p and Ilj. 

The procedure of the ISM for the initial value problem of the matrix NLS equation 
( 12. ip is summarized as follows. First, we solve the eigenvalue problem (13.11) for the initial 
value Q{x,0) and obtain the scattering data at t = (direct problem). Then, with the 
time-dependent scattering data (13.851) . we solve the Gel'fand-Levitan-Marchenko equa- 
tions (13.671) and (13.691) to obtain IC{x,y,t) and Q{x,t) (inverse problem). This procedure 
provides the direct proof of the complete integrability of the matrix NLS equation (12. ip 
under the nonvanishing boundary conditions (12. 6p . Further researches are required to 
establish each step in a rigorous way. 

3.4 Soliton solutions 

We construct soliton solutions of the matrix NLS equation under the reflection-free condi- 
tion: 

B{X,0=B{\,0=O, (e:real), (3.87) 
whereby the first terms, J-'c{z,t) parts, in eqs. (13.86P identically vanish, 

N N 

J^^{z, t) = -J2 i^i(t)e^^^^ ^2(2, t) = -Yl ^^j'^jity'-'"- (3-88) 
and similarly for and JF2. Assume the form, 

N 



ICu{x,y,t) = iAo5^/cSn,(t)e^^^(^+^), (3.89) 

N 

]Cu{x,y,t) = iAo5^/ci^n,(t)e^^^(^+^). (3.90) 
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Then, the Gerfand-Levitan-Marchenko equations fl3.67p and fl3.69p are reduced to a set 
of hnear algebraic equations, 



iAo^g + i(0 - A,)/ + Y: ^#^7^^- n,(t)e'c- + Y: 7T^^!?n,(t)e^c.. = o 

^ HO + CfcJ ^ HO + Oj 



(3.91) 



for j = 1, ■ ■ ■ , A^. Those are simphfied into 



'^11 
-/ 



0-A>-' 



(3.92) 



(3.93) 



A, 



N 



~ A,- 



^ i(o + Cfe)(a-Afe) 



:tKO + o)(o-A,) 



(3.94) 

Replace Aj — > — Aj for notational reason, and note the relations f l3.75p and that the time- 
dependence of Uj and flj is 

Ujit) = nj(o)e^% Ujit) = ^J■(o)e'^^ (3.95) 

where we introduce the function 

Xj{x, t) = 2iO(x - 2(Aj + k)t). (3.96) 

Then, from eq. (13.940 we have a matrix form, 

•51,1 ■ 



'^12 ■ ■ ■ '^12 



Sn,i 



Sl,N 

Sn,n 



N 

(iW), 



(3.97) 



where for j = 1, ■ ■ ■ , A^, 



Si/ 



-6, J- 



1 1 

+ 



n,;e> 



(3.98) 



+ A, i(0 + 0) V + Ai + A, . 
Solving the linear equations fl3.97p and using eq. fl3.49ap . we arrive at the multi-soliton 
solution for the modified version of the matrix NLS equation (12.120 . 



N 



Q{x,t) = Ao/ + 2iAo(/---/)S^' 



Hie^i 



(3.99) 
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Before concluding, we summarize about the parameters of the solution. For e = +1, the 



solution fl3.99p is the A^-soliton solution. (j 



A^) is a real constant such that 



— Ao < Xj < Aq. Cj = i(-^o ~ A|)^/^ is pure imaginary. Ilj is an / x / Hermitian matrix. For 
e = -1, the solution (K99h is the M(= A^/2)-soliton solution. Xj and (j = (A^ + A^)^/^ 
(j = I,-- - ,N) are complex constants satisfying X2k = Xli^_^ and C2fc = —Qk-i f'^^ ^ — 



N/2. Uj is an / X / matrix satisfying 11 



2fc-l 



n 



t 

2k- 



Although, in a strict sense, we should impose that detllj = for all j, we can relax 
this condition. The reason is that in the limiting case where two distinct Aj and Xj merge 
in the (A^ + l)-soliton solution, the expression eq. ( 13.99^ for the solution remains true with 
replacing formally Ilj — Ilj + Uj. Recall that det(nj + Ilj) 7^ in general, even when 
det Hi = and det Uj = 0. 

Finally, multiplying the carrier wave part e^^'^^"'^*) as eq. fl2.10p . the multi-soliton solu- 
tion for the original NLS equation (12. ip under the nonvanishing boundary conditions (12.60 . 
(12. 9p is obtained: 



TV 



Q{x, t) = Aoe'^'^"-"*) <^ / + 2i(/ ■ • ■ I)S-' 



IlAre 



XN 



(3.100) 



As mentioned in Sec. 12. ![ the multi-soliton solutions for general nonvanishing boundary 
conditions are easily obtained through the unitary transformations (12. 8p . 



4 Concluding Remarks 

In this paper, we have studied both the self-defocusing and the self-focusing matrix nonlin- 
ear Schrodinger (NLS) equations under nonvanishing boundary conditions. Introducing the 
Lax pair, we have made the inverse scattering method (ISM) analysis for the systems and 
shown that the initial value problem is solvable. From the Gel'fand-Levitan-Marchenko 
equation with the reflection-free condition, the multi-soliton solution is obtained explicitly. 
The conservation laws are given in the same manner as for the vanishing boundary condi- 
tions [13], leading to an infinite number of the conserved quantities. On the other hand, 
the Lax pair, the direct and inverse problems, and the multi-soliton solution are altered 
due to the boundary conditions. 
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The multicomponent systems with nonvanishing boundary conditions include several 
important models in physics, for example the spinor model for Bose-Einstein condensates 
with repulsive and antiferromagnetic interactions. The equation for the dynamics of F = 1 
spinor condensate falls into the / = 2 case of the matrix NLS equation. As expected from 
the applicability of the matrix NLS equation to the spinor bright solitons [T¥[ IT5] . it is 
interesting to analyze the spinor "dark" solitons based on the results obtained in this 
work. We report this issue in a separate paper [T9] . 
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